Backward Raman compression of x-rays in metals and warm dense matters 
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Experimentally observed decay rate of the long wavelength Langmuir wave in metals and dense 
plasmas is orders of magnitude larger than the prediction of the prevalent Landau damping theory. 
The discrepancy is explored, and the existence of a regime where the forward Raman scattering is 
stable and the backward Raman scattering is unstable is examined. The amplification of an x-ray 
pulse in this regime, via the backward Raman compression, is computationally demonstrated, and 
the optimal pulse duration and intensity is estimated. 
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Coherent intense x-rays would enable various potential 
applications P, Q, however, it is very difficult to com- 
press an x-ray pulse, or even in the UV regime. Some 
progress has been made in this direction, based on the 
recent advances in the areas of free electron laser Q 
and the inertial confinement fusion It is shown that 
the backward Raman scattering (BRS) [1, @]; where two 
light pulses and a Langmuir wave interact one another 
to exchange the energy, is one promising approach to 
create an ultra short light pulse (down to a few atto- 
seconds 0, Hi); a pulse gets compressed and intensified 
via this laser-plasma interaction. The BRS is successfully 
used to create intense pulses of the visible light frequen- 
cies d, [l3| • It is natural to consider the same technique 
for the x-ray compression 0, However, some physical 
processes in the regime where x-rays might be compress- 
ible are considerably different from those in the visible 
light regime. There are new physical processes to be con- 
sidered, such as the Fermi degeneracy and the electron 
quantum diffraction to mention a few. 

We consider one key aspect for the BRS in the x-ray 
compression regime, i.e., the damping of the Langmuir 
wave. In an ideal plasma, the decay rate of the Lang- 
muir wave increases rapidly as the wavelength decreases. 
It poses a concern for the BRS, as the plasmon from the 
forward Raman scattering (FRS), having a lower wave 
vector than the BRS, depletes the pump. There have 
been various attempts to suppress the FRS in the visi- 
ble light regime [l^. The situation is different in metals 
and the warm dense matters that we consider for the 
x-ray compression. The interaction of electrons via the 
inter-band transition (the Umklapp process) [l6l - [2^ be- 
comes the dominant plasmon decay process, in contrast 
to the Landau damping in an ideal plasma. As a conse- 
quence, the decay rate of the long wavelength plasmon 
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is much higher than the prediction of the prevalent di- 
electric function theory. In this paper, we show that a 
parameter regime where the BRS is unstable and the FRS 
is stable does exist in metals and warm dense matters, 
as the plasmon from the FRS strongly decays. Here, we 
consider metals in room temperature, and show that an 
x-ray pulse can be compressed in this regime. We esti- 
mate an optimal pulse duration and intensity of the seed 
and pump pulses, with which the inverse bremsstrahlung 
would be prevented. We use an one-dimensional (1-D) 
simulation to demonstrate that a pump pulse with a du- 
ration of hundreds of femto seconds could be compressed 
to a pulse of a few or even sub-femto seconds. 

The three wave interaction can be described by p3l |: 
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where Ai = eEi/rrieUJiC is the amplitude of the pump 
(j — 1) and the seed {i — 2) pulse, Ei is the electric field, 
TOg is the electron mass, uji is the pulse frequency, c is 
the speed of the light, Vi is the group velocity of light 
(plasmon) scaled by c, Vi is the inverse bremsstrahlung 
rate, c; = uy^^/2LUi, and ki is the wave number of the 
pulses; A^ = fie/ne is the Langmuir wave amplitude, 1/3 is 
the plasmon decay rate, and C3 = (cg)^/2w3, where lo^ — 
cjpc = Aiml/nie, and q is the Langmuir wave number. 
In the BRS (FRS), the energy conservation reads loi = 
L02 + Wpo and the momentum conservation reads g = fci -t- 
^2 [q — ki ~ k2) ■ The plasmon decay, 1^3, and the inverse 
bremsstrahlung, i^i and V2, are to be discussed later. 

The plasma density response given by the third equa- 
tion in Eq.(IT]) is derived from a fluid equation. However, 
it can be still used in a degenerate or a partially degener- 
ate plasma by the following reason. The plasma response 
is given as ne{q,uj) = a{q, UJ3) / {l + ATre"^ / q'^a{q, uj3))4)pod, 
where 0pod is the ponderniotive potential from the beat- 
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ing of the seed and pump, and e = 1 + {Aire^ / q^)a (a) 
is the dielectric function (susceptibihty) , as an example 
of which, in a degenrate plasma, there is the well-known 
Lindhard function in a zero temperature or a finite tem- 
perature 24). It is shown by Capjack 2M t^^-t the fluid 
equation above describes reasonablly well the full kinetic 
plasma responce {n^) for the classical dielectric function 
if 2^3 = Im(e) — iRe(e). Since the theory to predict he 
in a degenerate quantum plasma is not different from 
the one in classical plasmas, the same fluid equation can 
be used as long as an appropriate dielectric function has 
been calculated. Here, we choose the uj^ in such a way 
that Re(e) = 0, then the damping, v^, is purely real, 
which we will assume. 

The Landau damping for the Langmuir wave, accord- 
ing to the dielectric function theory, is given as 

2 = ilmKg,^)], (2) 

where 7 is the decay rate, e is the dielectric function sat- 
isfying Re[e] = 0, and q is the wave vector. The Landau 
damping rate in the free-electron plasma can be com- 
puted by using either the classical dielectric function or 
the degenerate Lindhard dielectric function. The damp- 
ing rates from the experimental measurements and differ- 
ent theories are compared in Fig. [U in a range of q. The 
damping rate from the Lindhard dielectric function [2^ 1 
is essentially zero in this range of g, which is not vis- 
ible in this scale. We note that DuBois' theory [26| . 
accounting for the dynamical correlation, predicts 7 is 
proportional to icilkp)'^ in the long wavelength limit (as 
q goes to zero), where kp is the Fermi wave vector. While 
all these theories predict the damping rate to vanish in 
the long wavelength limit in a way or another, there are 
experimental evidences that the rates are finite in this 
limit. For example, experimental measurements in the 
(1,0,0) direction of Al exhibit finite damping rates as 
shown in Fig. [T] The experimentally observed plasmon 
decay rates do not vary much for q < 0.5fci?, where the 
inter-band transition is the dominant decay mechanism. 
It rapidly increases for q > 0.5kp (not shown in Fig. [T] 
see Fig. 3 in Ref. [13); where the Landau damping be- 
comes dominant. We express the experimental decay rate 
of the plasmon for q < 0.5kp as 

i^siq) = ?7(g)wpe, (3) 

where 77(g) = r/o + dr] / dq"^ {q'^ / kp) . For a typical metal, 
0.02 < ??o < 0.2 and drj/dq"^ ^ arjo, where 2 < a < 10 [H]. 
This experimental damping rate is successfully modeled 
by Sturm [l^], who takes into account of the Umklapp 
process of electrons in the presence of the spatially peri- 
odic ions. 

The inverse bremsstrahlung (vi and 1^2) in metals is 
given as [soj . 

Mlu,E,) - 47m,Z, (4) 
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FIG. 1: The damping rates of plasma of Ue — 1.8 x 10^^/cc; 
the prediction from the classical dielectric function (solid line; 
Te — 20eV, chosen to be comparable to the Fermi energy) 
and from the DuBois' theory (dashed line; corresponding to 
Te — 0). The experimental measurements for Al are adapted 
from Fig. 3 in Ref. [13], which are originally from Ref. [27l ] 
(crosses), Ref. [2^ (circles), and Ref. ,29] (triangles). 

where Ki — hu!i/2'meVp < 1, vp is the Fermi energy. 
Si = eEi/rrieOJiVp, of = 26"^ Ef / mfiwf , and F{a) can 
be approximated as for a <2 and (2/7rQ;) log(2a)^ 
otherwise. In Eq. it is assumed that hoj > and 
Huj > Ep, where Tg and Ep is the electron temperature 
and the Fermi energy, respectively. For a < 2, the above 
equation can be simplified as 

Ui = 2.9 X 10^"^ ZiK;^^^-^ sec-^ . (5) 

When 1 / Ui is shorter than the pulse duration, the pump 
pulse would get heavily damped before having a chance 
to compress the seed pulse. Note that 1^1 is proportional 

to (ujpc/uJl)^. 

From Eqs. (H]) and ([3]), the BRS growth rate, assum- 
ing that the pump intensity is large enough to make it 
unstable, is roughly estimated as 

9B = — Wpo. (6) 

The larger gs is, the stronger the BRS is. It can be seen 
from Eqs. (O and © that wi/wpo » 1 is desirable, as this 
condition would suppress the inverse bremsstrahlung, 
while keeping the BRS growth rate high. In order to 
maximize the efficiency, we choose Wi/wpc so that the 
wave vector (wi +uj2)/c does not exceed 0.25kp, beyond 
which the linear Landau damping is important. In vari- 
ous metals, the desirable uji/cupc is 

10 < — < 30, (7) 

which corresponds to the photon energy in AL between 
150 eV and 450 eV. 
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We obtain the possible operating regime for the Raman 
compression using Eqs. ([1]), ([3]), © and ([7]). The stability 
condition for the FRS, obtained from Eq. ([1]), is given to 
be C2C3|Ai|^ < i'2^'3, where V2 and can be found in 
Eqs. ([3]) and dU, and C3 — ujpd'^- This condition can be 
written as \Ai\ > \Aif\, where 



(8) 



and qp — Wpc/c. The same condition for the BRS is 
estimated to be \Ai\ < \Aib\, using C3 = {oji+oj2)'^ /2LUpc, 
where 



1^ 
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and qb = (wi + W2)/c. l^iBp is smaller than |^ifP by 
a factor of {oJi+oj2)'^ /ojpdr]{qp)/r]{qB)), which is consid- 
erably large since '>l{qF)/v{QB) — 1- 

A regime, where the BRS is unstable and the FRS is 
stable, exists due to the strong inverse bremsstrahlung 
and a rather constant damping of the Langmuir wave in 
a range of q. While these boundary is estimated assum- 
ing the maximum plasma wave amplitude using the large 
damping limit, this is useful for an estimation. Further- 
more, in the linear regime where the pump is uniform and 
the seed pulse is small so that the linear theory is valid, 
the above boundary is exact. Even when the FRS is un- 
stable, it would not be detrimental unless the background 
noise plasmon grows to such an extent that depletes the 
pump. The growth rate of the FRS, assuming unstable, 
is given as 



3^ = 0.25x1^^^1^1 1 

77 UJ2 



pe. 
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Unless gpTi > 10 , where ri is the duration of the pump, 
the FRS instability can be ignored. In AL, 77 is roughly 
0.03. Even for WpoTi = 10'^, the pump pulse intensity has 
to be unrealistically large for the FRS depletion to occur. 
In other words, even though the FRS is absolutely un- 
stable (linearly unstable), the FRS/RBS is convectively 
stable. 

Now a detailed estimation using the parameter for AL 
is in order, assuming Zi — 3, rj — 0.03, hujpc — 15.3eV, 
and Ep = 11.7eV, with which the electron density rig 
is 1.8 X lO^^/cc and 2.4 x lO^^sec"^ Consider 

the case when uji/ujpc ~ 10 (and 20). The inverse 
bremsstrahlung for the pump and the seed is roughly 
given as 1^1,2 = 7 x 10~^ti;pc (2^1,2 = 5 x 10~^a;po). The 
threshold for the FRS in Eq. ([5]) is estimated to be 
\Aif\'^ = 6.72 X 10--^ (l^iFp 1.19 X 10"*), which 
corresponds to 1.27 x lO^^ W/cm^ (8.98 x lO^^W/cm^). 
On the other hand, the estimated threshold for the BRS 
in Eq. © is l^isl^ = 1.68 x lO"'^, or 3.1 x lO^*^ W/cm^ 
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3.3 X 10^^ or 6.1 x lO^^^/cm^ 



We numerically integrate the 1-D model, Eq. ([T|), to 
demonstrate the compression of an x-ray pulse in AL 



(Fig. [2]). The pump pulse travels from the left to the 
right as does the seed pulse. The initial intensity of the 
pump is chosen to be 2.5 x 10^^ W/cm^ in Fig. [5] (a) and 
10^8 W/cm2 for Fig.[2](b). The seed is lO^^ W/cm^ and 
the pump frequency is wi/wpe = 15 in both cases. The 
seed pulse extracts the energy from the pump and the in- 
tensity gets amplified by a factor of about 15 (Fig. [5] (a)) 
and 140 (Fig.[5](b)). When the nonlinear pump depletion 
is negligible, tb = l/gB m Eq. (jB)) is the growth time of 
the seed pulse, which roughly corresponds to the width of 
the compressed pulse. When the nonlinear pump deple- 
tion becomes important, the pump depletion time might 
be shorter than the linear growth time tb by a factor 



of |^2Mi 



Since the pump feeds the energy into a 



small portion of the seed pulse in the nonlinear deple- 
tion regime, the actual compressed pulse width would be 
much narrower than the prediction of our linear analysis 
in Eq. In Fig. [2] (a), the maximum of the Langmuir 
wave intensity {A^l during the compression remains less 
than about 0.3. 

For an intense Langmuir wave, the nonlinear cascading 
could be an important saturation factor (SlJ. It is esti- 
mated that the nonlinear cascading is not important in 
the compression time scale (tc = lO/wpo). In addition, 
the damping of the Langmuir wave is strong as shown 
here, and the damping of the ion acoustic wave is strong 
due to the Umklapp process. Consequently, the mini- 
mum threshold for the nonlinear cascading [31| would be 
high, so it becomes of less concern in dense plasmas. 

Now we consider the wave breaking. The maximum 
strength of the Langmuir wave is |^3|max = |^i||^2|(wi-l- 
uj2y /i^pcVil) (the third equation in Eq. ([T])). The theo- 
retically achievable value of |^3|max is less than the unity, 
which is called the wave breaking. Before 1^3 1 reaches the 
wave breaking, various saturation mechanisms come into 
play, including the frequency detuning and the soliton 
formation [3^ . [ssj . These mechanisms make the linear 
theory presented here become insufficient. In the simula- 
tion results shown in Fig. [21 the Langmuir wave frequency 
is higher than Ep/h and k = 0.2kp, and the phase veloc- 
ity of the plasmon is much higher than the Fermi velocity. 
We estimate that, until IA3I reaches ~ 0.5, the frequency 
detuning due to the trapped particles is negligible and the 
above linear theory is valid, which is the case in Fig. [5] 
(a). When I A3 1 max exceeds ~ 0.5, the currently available 
theories [32 1 need to be extended by taking into account 
of the electron deg eneracy and the quantum diffraction 
effects [nl, [H [ij In Fig. H (b), we assume that the 
1^3! quickly decays when {A^l > 0.5 (due to the wave 
breaking, the frequency detuning, and the nonlinear cas- 
cading) . 

The inverse bremsstrahlung would be of less concern 
as can be seen from Eq. and it would be further 
diminished by the electron degeneracy [llj, ■ As the 
pump travels inside the medium, it heats electrons via 
inverse bremsstrahlung especially on the back of the 
pump. As the electron temperature increases, the in- 
verse bremsstrahlung decreases faster. While a high elec- 
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FIG. 2: Time series of the seed pulse, A2 = eE2/me(jJ2C, and 
the Langmuir wave, A3 = h/n, traveling to the right through 
AL (gray area, < a; < 2000), where x is normalized by ujpa/c. 
Each pulse is separated by SOOO/upe, and the duration of the 
pump pulse is 4000/ajpo in both cases. See the text for details. 



tron temperature would certainly increase the Landau 
damping and the partial wave breaking arising from the 
trapped particles, it is suggested that the damping and 



the wave breaking would be weak in dense plasmas fl^ . 
More concern for the heating of electrons would be the 
detuning effect of the plasmon from the change of the 
plasmon disperson relation, and the reduction of the plas- 
mon damping. The reduction of the plasmon damping 
would make the FRS more unstable, but our preliminary 
estimation suggests that this is still containable. The 
more detailed estimation should wait until the thoory of 
the plasmon damping for the partially degenerate plasma 
is available, which is in progress. 

To summarize, it is shown that the FRS is less harm- 
ful for the Raman compression in metals and warm dense 
matters than in the gas discharge plasmas or ideal plas- 
mas, as a consequence of the enhanced decay of the long 
wavelength Langmuir waves. The optimal frequency and 
the pulse duration for the compression is estimated in our 
analysis. Our 1-D simulation based on the linear anal- 
ysis demonstrates that a pump pulse with the duration 
of lO'^/wpo — lO'^/wpo and the intensity comparable to 
lO^^W/cm'^ is compressed to the one of l/wpc — 10/wpo, 
with the intensity of lO^^W/cm^ - lO^iW/cm^. 

As suggested above, it is crucial to understand the 
plasmon damping rate in the presence of the relic lattice 
structure. The plasmon decay in metals strongly depends 
on the direction relative to the lattice structure. As the 
electron temperature increases, the electron screening of 
the ions would get diminished, which enhances the plas- 
mon damping [T3|. As the electron heats up further, the 
ions would finally lose their lattice structure. A larger 
value of 77(g) in Eq. ^ would be preferred for a strong 
pump, while a small value of 7y(g) is better for a weak 
pump. The damping rate in the warm dense matters 
might be accessible by the measurement of the electron 
energy loss in the thin heated foil experiments [l^. The 
effect of the phase transition on the plasmon decay is the- 
oretically challenging to understand, yet it is important 
because of its relevance to the BRS x-ray compression. 

The authors would like thank Prof. Fisch for useful 
discussion and advice. 



[1] H. N. Chapman, Nature Physics 2, 839 (2006). 

[2] M. Hentschel, R. Kienberger, C. Spielmann, G. A. Rei- 

der, N. Milosevic, T. Brabec, and P. Corkum, Nature 

414, 509 (2001). 
[3] P. Emma, K. Bane, M. Cornacchia, Z. Huang, H. Schlarb, 

G. Stupakov, and D. Walz, Phys. Rev. Lett. 92, 074801 

(2004). 

[4] M. Tabak, J. Hammer, M. E. Glinsky, W. L. Kruerand, 

S. C. Wilks, J. Woodworth, E. M. Campbell, and M. J. 

Perry, Physics of Plasmas 1, 1626 (1994). 
[5] V. M. Malkin, G. Shvets, and N. J. Fisch, 

Phys. Rev. Lett. 82, 4448 (1999). 
[6] N. J. Fisch and V. M. Malkin, Phys. Rev. Lett. 82, 4448 

(1999). 

[7] V. M. Malkin and N. J. Fisch, Phys. Rev. Lett. 99, 
205001 (2007). 



[8] V. M. Malkin, N. J. Fisch, and J. S. Wurtele, 

Phys. Rev. E 75, 026404 (2007). 
[9] Y. Ping, W. Cheng, S. Suckewer, D. S. Clark, N. J. Fisch, 

M. S. Hur, and J. S. Wurtele, Phys. Rev. Lett. 92, 175007 

(2005). 

[10] W. Cheng, Y. Avitzour, Y. Ping, S. Suckewer, N. J. 

Fisch, M. S. Hur, and J. S. Wurtele, Phys. Rev. Lett. 

94, 045003 (2005). 
[11] S. Son and N. J. Fisch, Phys. Rev. Lett. 95, 225002 

(2005). 

[12] S. Son and N. J. Fisch, Phys. Lett. A 329, 16 (2004). 
[13] S. Son and S. Ku, Phys. Plasmas 17, 024501 (2010). 
[14] S. Son and S. Ku, Phys. Plasmas 17, 010703 (2009). 
[15] V. M. Malkin, Y. Tsidulko, and N. J. Fisch, 

Phys. Rev. Lett. 85, 4068 (2000). 
[16] K. Sturm, Z. Physik B 25, 247 (1976). 



5 



[17] K. Sturm, Z. Physik B 28, 1 (1977). 

[18] P. C. Gibbons, S. E. Schnattcrly, J. J. Ritsko, and J. R. 

Fields, Phys. Rev. B 13, 2451 (1976). 
[19] W. Ku and A. G. Eguiluz, Phys. Rev. Lett. 82, 2350 

(1999) . 

[20] W. Ku and A. G. Eguiluz, J. Phys. Chem. Sohds 61, 383 

(2000) . 

[21] J. M. Pitarke and I. Campillo, J. Phys. Chem. Solids 164, 

147 (2000). 

[22] G. Tirao, G. Stutz, V. M. Silkin, E. V. Chulkov, and 
C. Cusatis, J. phys.: Condens. Matter 19, 046207 (2007). 

[23] C. J. McKinstrie and A. Simon, Phys. Fluids 29, 1959 
(1986). 

[24] J. Lindhard, K. Dan. Vidensk. Sels. Mat. Fys. Medd 28, 
8 (1954). 

[25] C. E. Capjack, C. R. James, and J. N. . McMuUin, 



,]. Appl. Phys. 53, 4046 (1982). 
[26] D. F. DuBois and M. C. Kivclson, Phys. Rev. 186, 409 
(1969). 

[27] P. J. Zacharias, J. Phys. F: Metal Phys. 5, 645 (1975). 
[28] M. Urner-Wille and H. Raether, Phys. Lett. 58A, 265 
(1976). 

[29] T. Kloos, Z. Physik 265, 225 (1973). 
[30] S. Son and N. J. Fisch, Phys. Lett. A 356, 72 (2006). 
[31] R. P. Drake and S. H. Batha, Phys. Fluids B 3, 2936 
(1991). 

[32] H. Rose and D. A. Russel, Physics of Plasmas 8, 4784 
(2001). 

[33] G. J. Morales and T. M. O'Neil, Phys. Rev. Lett. 28, 417 
(1972). 



